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and find that such representations exist in terms of operators that are state- 
dependent in a characteristic manner. 
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I. INTRODUCTION 



The selection of collective variables is an important problem in the study of large ampli- 
tude collective motion. In the constrained Hartree-Fock (CHF) or Hartree-Fock-Bogoliubov 
(CHFB) calculations, the collective subspaces are generated by a small number of one-body 
constraint (also called generalised cranking) operators which are most commonly taken to 
be multipoles that can be represented as homogeneous polynomials in the coordinates, i.e., 
t'^Ylk- In realistic calculations of processes such as fission, the number of coordinates needed 
to describe the full nuclear dynamics can easily become larger than what can be dealt with 
in satisfactory manner, and a method to determine the optimal combination needs to be 
devised. Furthermore, there is no a priori reason to limit oneself to multipole operators, 
and the cranking operators should be determined by the nuclear collective dynamics itself, 
from amongst all possible one-body operators. 

In our past work, we have investigated a theory of adiabatic large amplitude collective 
motion as a method to generate such self-consistent collective subspaces (see Ref. and 
references therein). The key ingredient of the method is the self-consistent determination 
of the constraint operator, and as such it may provide an answer to the selection question 
discussed above. Using the local harmonic version of the theory the local harmonic 
approximation (LHA), we have recently embarked on a study of the properties of large 
amplitude collective motion in systems with pairing. So far we have dealt with two simple 
models: The first one is a semi-microscopic model of nucleons interacting through a pairing 
force, coupled to a single harmonic variable 0. The second is a fully microscopic 0(4) 
model which may be regarded as a simplified version of the pairing-plus-quadrupole (P+Q) 
Hamiltonian studied in this paper. It has turned out that the self-consistent collective 
coordinate obtained by the LHA accounts quite well for the exact dynamics of the models. 
We have also shown that the CHFB calculations using the mass-quadrupole operator as the 
constraint operator can result in incorrect results 0. 

In this paper, which is an attempt to move towards realistic nuclear problems, we will 
encounter some problems that need to be resolved first: The main problem is the treat- 
ment of conservation laws. In studies of large amplitude collective motion, the mean-field 
states often do not have the symmetries of the Hamiltonian. Then, the spurious (Nambu- 
Goldstone) degrees of freedom, such as the translation and rotation of nucleus, inevitably 
appear, which must be properly separated out from the other collective degrees of freedom 
like nuclear shape change. In the models discussed in Refs. |@,|^, the spurious degrees of 
freedom could be removed by hand thanks to the simplicity of the models, however in re- 
alistic nuclear problems, it is almost impossible to do so. The second and more obvious 
problem is the increase in the number of degrees of freedom. Within the time- dependent 
Hartree-Bogoliubov (TDHB) approximation, the dynamics of the models of Refs. is 
described by only a modest number of degrees of freedom (4 ~ 12). For realistic problems 
in heavy nuclei, on the other hand, we need to deal with at least tens of thousands, and 
possibly millions of degrees of freedom! These two points, the effect of the conservation laws 
and the applications to realistic problems, are the main issues discussed in Sec. y. 

In the LHA, the collective path (or collective manifold for more than one coordinate) is 
determined by solving the CHFB problem with a cranking operator which is self-consistently 
determined by the local random-phase approximation (RPA). It is of course a well-known 
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fact that the spurious modes decouple from the physical modes in the RPA. However this 
is true only when the system is at a minimum of the mean-field potential. Since the local 
RPA needs to be solved at non-equilibrium points as well, we need to study the properties 
of the LHA equations. It will be shown that the symplectic version of the LHA combined 
with an extended adiabatic approximation properly takes account of conservation laws. It 
will also be shown that this formalism naturally explains the role of the chemical potential 
term — AA^ in the Hamiltonian of the quasiparticle RPA. 

Since the LHA procedure requires one to solve the local RPA at large number of points 
along the collective path (manifold), the large dimension of the RPA matrix for realistic 
problems leads to very heavy computational burdens. Thus, instead of solving the local 
RPA, we would like to preselect a small group of operators, and find the RPA modes as an 
optimal combination of these operators at each point of the collective surface. This means 
that we restrict the RPA diagonalisation to this small space, rather than deal with the full 
millions-by-millions RPA matrix [Q. In this paper, we attempt to find a set of one-body 
operators which approximate well the self-consistent cranking operator. The pairing-plus- 
quadrupole (P-l-Q) Hamiltonian is adopted for the analysis because the RPA calculation is 
relatively easy due to the separability of the force (some of the results for the Sm isotopes 
have already been reported in Ref. [Q). Since the P+Q model is known to be able to describe 
collective phenomena involving both pairing and quadrupole degrees of freedom we 
expect that a similar choice for the set of one-body operators should work for other realistic 
Hamiltonians as well. 

The P-l-Q model is probably one of the most simple and successful nuclear Hamiltonians 
which allows us to discuss realistic problems involving pairing and quadrupole degrees of 
freedom. Baranger and Kumar analysed in great detail the (adiabatic) collective motion in 
the P-l-Q model assuming that the collective variables are the mass quadrupole operators 
||T0HT4| . Thus, they reduced the large number of two-quasiparticle (2qp) degrees of freedom 
(for this model of the order of a thousand) to only two "collective" coordinates, (3 and 
7. However, our previous study of the 0(4) model suggests that even for such simple 
Hamiltonians, the self-consistent collective coordinates are often not as simple as one might 
think. The P-l-Q model (with an additional quadrupole pairing interaction) has also been 



studied by Kishimoto and Tamura ||T5|,|T6[ using boson expansion techniques. They found 
that the coupling to (non-collective) 2qp excitations is very important, which also suggests 
that the collective coordinates are not just the mass-quadrupole ones. In this paper we 
show that the normal-mode coordinate of the random-phase approximation is quite different 
from the mass-quadrupole operator, and how the structure of the self-consistent collective 
coordinates changes as the nucleus is deformed. 

In Sec. H, we give an account of the general theoretical arguments. The problem of the 
spurious modes and the extension of the RPA formalism are discussed in Sec. |II B| and |II C 



In Sec. [II m , the method of projection of the RPA equation onto a set of one-body operators 
is described. Numerical calculations for the RPA and the projected RPA for the P-l-Q model 
are given in Sec. |T| and finally we give the summary and outlook in Sec. |V|. 
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II. LOCAL HARMONIC FORMULATIONS AND CONSERVATION LAWS 



In this section, we use a summation convention where the repeated appearance of the 
same symbol for upper and lower index indicates a sum over this symbol for all allowed 
values. We also use the convention that a comma in a lower index indicates the derivative 
with respect to the coordinate, i.e., = OF/ 8^°". 



A. The TDHFB theory and classical constants of motion 

Utilising the TDHFB theory we transcribe the original quantum Hamiltonian into the 
classical Hamiltonian. We follow Ref. |T7] to define the canonical variables. We denote 
a time-dependent generalised Slater determinant by |^l'(t)) which is characterised by the 
(quasi-) density and pair matrices. 



p,j = (^(t)|a]ai|^(t)), = {^{t)\aja^\^{t)), 



(2.1) 



where |\E'(t)) is assumed to be normalised at all times ((\Ef(t)|\E'(t)) = 1) and the quasiparticle 
operators (aj,a|) are defined with respect to a reference state I^E'o)- Based on a classical 



Holstein-Primakoff mapping ]TB| , p!U 



Pij 



(2.2) 



the density and pair matrices are now mapped onto complex canonical variables Pij. As 
usual real canonical variables ^ and vr are introduced as P = + in) / \/2. In this canonical 
parametrisation, the reference state l^'o) corresponds to the origin of the phase space, where 
^ = vr = 0. 

We shall specifically look at one-body hermitian operators R and anti-hermitian opera- 
tors 



J2S'{i3)alaj, 



(2.3) 
(2.4) 



i>j 



where Rq = (\&o|-R|^o) and all the matrix elements, R{ij), R'{ij), S{ij) and S'{ij), are 
assumed to be real. Here R'{ij) {S'{ij)) are symmetric (anti-symmetric) with respect to 
permutation of the indices i and j. The corresponding classical representations, 71{C,, tt) and 
iS(^,7r), can be found by replacing the fermion pair operators a]aj and ajOj by the density 
and pair matrices pij and Rij of Eq. ( p.2|) . Expanding p and R with respect to vr, we find 
that TZ{^, tt) consists of even-order terms in tt while S{^, vr) contains odd-order terms only. 



7^(e,7^) = 7^W(0 + -7^(2)-^(Ov^„v^^ + 0(7r^), 
5(e,7r)=5«"(07r, + O(7r3). 



(2.5) 
(2.6) 



Starting from these equations, pairs of two-quasiparticle (2qp) indices {ij, kl, ■ ■ ■) are denoted 
collectively by a single Greek index (a, (3, ■ ■ ■). The 2qp matrix elements with respect to 
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the reference state, R{ij) and S{ij) which are now denoted as R{a) and S{a) respectively, 
are related to the derivatives of the classical representations. 



R[a) = —j= 



V2 di- 

i dS 
b[a) = —j= 



C=7r=0 



\f2 dTXa 

The classical Hamiltonian is given by 



V2 



(2.7) 
(2.8) 

(2.9) 



where the right hand side is calculated in terms of the density and pair matrices ( |2.2| ). The 
expansion of 7i(^, vr) in powers of vr defines the potential V{^) = Tiiji = 0) (zeroth order) 
and the mass parameter B°'^{S,) (second order). 



5"^(0 



(2.10) 
(2.11) 



7r=0 



The terms of order of tt^ and higher in the Hamiltonian are neglected in the adiabatic 
time-dependent mean- field theory. The tensor -Bq,/3, which is defined as the inverse of 
[B°''^B^I3 = 6"^), plays the role of metric tensor in the Riemannian formulation of the local 



harmonic approximation (LHA) as discussed in Sec. ilU . 

It is well-known that the TDHFB equation preserves the symmetries of the quantum 
Hamiltonian: If the one-body operator P commutes with the Hamiltonian, [P, H] = 0, the 
classical representation of P, V{^,7i) = (\E'(t)|P|\l/(t)), is a classical constant of motion. 
Consequently, the Poisson bracket between V and 7i must vanish. 



{v,n} 



PB 



0, 



which can be used to show that 



-p(O) ^aP _ p(2)al3y 



0, 



(2.12) 



(2.13) 
(2.14) 



which are the terms of zeroth and first order in vr in Eq. ( pJ2[) , respectively. Here, 



and p(2)«/3 

are defined, according to Eqs. ( p.5|) and ( p. 61) , as terms of the zeroth, first 
and second order in vr, respectively. The equations ( |2.13| ) and ( |2.14| ) hold at arbitrary points 
in configuration space. 



B. Riemannian version of the local harmonic approximation (LHA) 

The adiabatic approach to large amplitude collective motion is based on a search for the 
collective (and non-collective) coordinates by performing a point transformation of the 
original coordinates 
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q' = nO, C = g''{q) (/i,« = l,---,n). (2.15) 
The conjugate momenta are given by 

Pi^ = = f%Pt^- (2-16) 
The adiabatic Hamiltonian ( p.lO| ) is transformed into 



n = V{q) + ]^B^''p,p, + 0{p'), (2.17) 
V{q) = V{giq)), B^^ = f^B^^f^'^. (2.18) 

The point transformation on the adiabatic Hamiltonian gives the same result as the one 
obtained by first performing this transformation on the exact Hamiltonian and then taking 
the adiabatic limit. This is not true for a general canonical transformation, specifically for 
the extended adiabatic transformation to be discussed in the next section. 

We divide the new coordinate set into three subsets, the collective coordinates q^, 
the spurious coordinates g^, J = 2, • • • , M + 1, and the non- collective coordinates g", a = 
M + 2, ■ ■ ■ ,n. For simplicity, we assume here a single collective coordinate. The spurious 
coordinates are related to the symmetry breaking in the mean field (Nambu-Goldstone 
modes). The local harmonic approximation consists of two sets of equations, the force 
condition and the local RPA, which define the collective manifold in configuration space (see 
Ref. for a complete description). When spurious modes are present, the force condition 
is modified slightly, 

K« = Ai/i + A//,!, (2.19) 

where the second term in the right-hand side is the additional term. The local RPA equa- 
tion can be formulated in different ways [H, and in this section we discuss the case of the 
Riemannian formulation, based on 

V^,B''^f}p = mffl, (2.20) 
glV.c.^B'^^ = mfg^- (2.21) 

Here the covariant derivative is defined by 

V,^p = V^p - Tl^V,, (2.22) 

where the affine connection F is defined with the help of metric tensor Baf^ as 

r)37 = 7^B"^ (^.5/3,7 + Bs-y,f3 - Bp^i^s) . (2.23) 



The collective path will be determined by solving Eqs. ( p.l9| ) and (|2.20|) selfconsistently. It 



consists in finding a series of points where the local RPA eigenvector f \ satisfies the force 
condition at the same time. 

The question at present is whether the solution of Eqs. ( ^.2(j| ) and ( |2.21D is orthogonal to 



the spurious modes. The spurious coordinates are normally linked to the one-body operators 
Pi which commute with the Hamiltonian. It is convenient to discuss separately the case in 
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which the symmetry operators Pj are hermitian with real matrix elements and the case in 
which Pj are anti-hermitian with real matrix elements. First, let us discuss the latter case. 
As is seen in Eq. ( p.6| ), neglecting the higher-order terms in vr, the classical representation 
of the symmetry operator is a classical momentum variable pj that is a constant of motion, 



p, = (vl/|P,|vI/)=p«%„ + 0(vr^). 



(2.24) 



From this, we can immediately see that g'} = vf^°' . Differentiating Eq. (|2.13|) with respect 
to keeping in mind that 'p(^)" must be proportional to (7°, we find 



Using an identity for the coordinate transformation ||T| 



+ J- 1319,119, u J- 



fii/ 



9:x 



this leads to 



9jVa,B^^ 



(2.25) 



(2.26) 



(2.27) 



where T is defined in the same way as F, ( |2.23| ), with obvious change of B to B. At the 
equilibrium = = 0, and since the right-hand side of Eq. (|2.27| ) vanishes, the spurious 
modes are indeed the zero-energy solutions of the RPA equation. Thus, all finite-energy 
solutions are automatically orthogonal to the spurious modes. However, at non-equilibrium 
points, the (7" are no longer RPA eigenmodes and there is no guarantee that the collective 
mode g"^ obtained from the local RPA equation ( |2.21| ) is free from spurious admixture]^. In 
Ref. 0, the local RPA equation was modified to insure that = = 0. However, in 

this formulation, we need to calculate both and (yf", which means that we have to solve 



the Thouless-Valatin equations for the spurious modes |22,23 



In the second case, where the symmetry operator P/ is hermitian, we have coordinates 
^ which are the constants of motion: 



q' = {%Pi\^) 



(2.28) 



If we limit ourselves to point transformations, we find Vf\Cj = /^(O ^'^d should ignore 
all terms depending on vr. At equilibrium Eq. ( |2.14| ) implies that the spurious modes 
correspond to zero eigenmode of the mass matrix P"^/^ = 0. Therefore we have a zero- 
energy mode, 



(2.29) 



Away from equilibrium, however, the are no longer RPA eigenmodes, due to the second 
term in (|2.14|) which is related to the terms we have neglected in Eq. (|2.28| ). 



1 If B^^^ and F^^, are block-diagonal for the spurious and the collective spaces on the collective 
path, the spurious solutions correspond to finite energy eigenvalues but are orthogonal to the 
collective mode. This is the case for HF problem in -^^Si as discussed in Refs. ij2|,|2|. 
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In summary, the Riemannian LHA can automatically separate out the spurious modes 
only at equilibrium = 0. Since the local RPA at equilibrium is nothing but the con- 
ventional RPA, the decoupling of the spurious modes is well-known. However, at non- 
equilibrium points, the spurious modes can in general mix with the collective one. This 
undesirable mixing disappears in the symplectic LHA with an extended adiabatic treatment 
described in the next section. 



C. Symplectic LHA and extended adiabatic transformation 

The symplectic version of the LHA is formulated with an affine connection 

f = ^?"m/S7' (2-30) 
which can be written in the usual form ( |^.23| ) by replacing the metric tensor by 

K^p^Y.f%f'p- (2-31) 

Since this metric depends on the final coordinates q^, we cannot define it from the beginning 
in contrast to the Riemannian metric Bap. In this formulation, the configuration space is 
assumed to be fiat and diagonal {K^j_i, = 6^u) in the final coordinate system q^. The covariant 
derivative is defined by 

= Vap - f%pV,- (2.32) 

For the case that the momenta pi are constants of motion, we can prove that the spu- 
rious modes are the zero-energy solutions of the symplectic local RPA equation. Indeed, 
differentiating the chain relation g^^fj^ = 5^ with respect to q^ , we obtain 

glJ'p = -Cglf'pr (2.33) 

Utihsing this equation and the definition of the covariant derivatives ( p.32|) , Eq. (|2.25| ) can 
be rewritten as 

= 0. (2.34) 

This is a consequence of the TDHFB equation of motion. Therefore the spurious modes are 
zero-energy solutions of the symplectic LHA equations, anywhere in the configuration space. 

In the case that the coordinates q^ are constants of motion, we need to lift the restriction 
to point transformations, based on the power expansion with respect to vr keeping the 
adiabatic assumption. Instead of Eqs. ( |2.15| ) and (|2.16| ), we generalise these to 

= no + ^/^'^""^(Ovr.vr^ + 0{n'), (2.35) 
e = g'^iq) + ^^^'^"^"^(^Kp^ + C?(/), (2.36) 

and 
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p^ = ^7>. + 0K), (2.37) 

7r. = />^ + 0(p=^), (2.38) 

where the terms cubic in momenta do not play a role in the modification of the theory. The 
original Hamiltonian (|2.1CI|) is transformed to, up to second order in p only, 

H{q,p) = Viq) + h^''p,p,, (2.39) 

= f^^B'^'^f'^^ + (2.40) 

Substituting Eq. ( |2.35| ) in Eq. ( p. 361 ), using Eq. ( p. 38] ), we find the relation 



g'-'^'^''' 9^91 = -f'^^^'9% (2.41) 
From the canonicity condition {q^, (f}'PB = 0, we also find 

The major difference between the use of the extended adiabatic transformation and a 
point transformation is the modification of mass parameter, 

= B"^ - (2.43) 



Here we have used the relation (p.41|) to obtain the last equation. The local RPA equations 



have the same forms as the Riemannian RPA equations, (|2.2CI|) and ( p.21|) , after replacing 
Vaf3 by Vap and 5"^ by 5°^. 

Let us return to the spurious modes. For hermitian constants of motion, i.e., coordinates 
q^ as in Eq. (|;2|), we can identify f{0 = ^f^(0 and f^^^^""^ = V^"^ for the extended 



adiabatic transformation. Utilising the TDHFB equation (|2.14|) and the canonicity condition 
(g), we find 



= B-^f;^ - f^'^'-^VJ^p 

= 0. (2.44) 

Therefore, anywhere in the space, it is guaranteed that the spurious modes are zero-energy 
solutions of the symplectic LHA with the extended adiabatic transformation. 

It is illustrative to discuss the difference between the Hartree-Fock (HF) and Hartree- 
Fock-Bogoliubov (HFB) approaches. The HF state is at an equilibrium point of the potential, 
Va = 0, where the force condition ( ^.19[ ) is satisfied with Ai = A/ = 0. We also see 



V;af3 = V;ai3 = V,ai3 and 5"^ = B"^ . Thus, any version of the LHA is equivalent to the 
conventional RPA. On the other hand, the HFB state is not a real equilibrium. We have 
a spurious coordinate associated with the particle number q^ = (\i/|A^|\i/) and the gradient 
of the potential has a non-zero component along the direction of q^ , Vn 7^ (the other 
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components of the gradient are all zero). The force condition is again satisfied with Ai = 
and Xn = Vn, but neither V^p = V^p nor = hold any more. Instead, we have 

^»I3 = 5"/3 _ Y^f{l)Naf3^ ^2.45) 

Vap = - ^iv/J^. (2.46) 

Since Vn is nothing but the chemical potential A a? at the (constrained) HFB minimum, 
the symplectic version of the LHA with the extended adiabatic approximation is equivalent 
to the conventional quasiparticle RPA with the constrained Hamiltonian, H' = H — X^N. 
This guarantees the separation of zero-energy spurious modes, while this is not the case for 
the Riemannian version of the LHA. This symplectic formulation also illuminates why the 
quasiparticle RPA at the HFB state should be performed using the Hamiltonian H', not the 
original H. 

As we have seen above, the symplectic formulation of LHA with the extended adiabatic 
transformation has an advantage over its Riemannian version. Unfortunately we do not know 
a general method to calculate /(^)/^"'^. It may be possible to calculate this quantity if the 
collective coordinate is explicitly given by a combination of elementary one-body operators 
(see the discussion in Sec. [rV| ). 

D. Projected LHA and the spurious modes 

One of the main obstacles in applying the LHA techniques described in previous sections 
to realistic nuclear problems is the fact that we need to diagonalise an RPA matrix of large 
dimensionality at each point of the collective path. In this section, we describe a method 
to approximate the RPA eigenvector by taking linear combinations of preselected one-body 
operators, generalising the method discussed in Ref. 0. 

First we select a small number of one-body operators, k = 1, ■ ■ ■ ,n' , assuming that 
the collective coordinate can be approximated locally as a linear combination of these 
operators, 

n' 

q'^i^im, f = J2ckF^'\ (2.47) 

k=l 

This means that the RPA eigenmodes f^a is projected onto the subspace spanned by {^q-*}, 

n' 

U^la = T.^kJ'^'\ (2-48) 

k=l 

where J-'^''^ are the classical representations (expectation values) of F^''\ In order to deter- 
mine the coefficients Ck in the linear combinations, the local RPA equation is also projected 
onto the subspace {^q-*}: 

M'^^ci = {hQfN''^ci, (2.49) 
where kCl is a frequency of the projected RPA and 
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y^kl ^ jr{k) QaPy^^^Q^S jril) ^ (2.50) 
^kl ^ jr(k)-ga(3jrW^ (2.51) 

Normally, in order to obtain the RPA eigenmodes and frequencies, we need to diagonalise 
the RPA matrix B°''^V^i3 whose dimension is equal to the number of active 2qp degrees of 
freedom. The dimension of matrices M'^' and N'^' is equal to the number of selected one- 
body operators {F^''^}. Therefore, if we can approximate the RPA eigenvectors by using a 
small number of operators, it will significantly reduce the computational task. 
The equation for the force condition is also projected the same way, 

yfe = AiN'^^Q, (2.52) 

with 

yk ^ jr(k)^a/3y^_ (2.53) 



Equations ( p.49|) and (|2.52|) must be solved self-consistently in order to determine a collective 
path. 

In cases where spurious modes exist, the spurious components should be removed from 
each operator as 

= - (2.54) 

and the force condition ( p.52| ) should be also modified by adding the terms A/W'^^ in the 
right-hand side with W'^^ = /"(^^^E^^/jj. 

The removal of spurious modes in Eq. (|2.54|) is a tedious task because we need to obtain 



both and g'} by solving the Thouless-Valatin equation. However, in case of hermitian 
constants of motion , adopting the extended adiabatic approximation in the previous 
section, we show that the explicit removal as in Eq. ( p.54|) is unnecessary. Using the modified 



mass parameter ( |2.43 ), we have B'^^f^p = 0, Eq. ( ^.44 ) showing that the correspond to 



the zero eigenmodes of mass matrix ( p.44|) . Therefore, we find 

M^' = J^(„^)5°^K^^5^^J^y = P^^B'^'^Vp^B^^^f, (2.55) 
j^ki ^ jr(fc)5"/3j-(0 = /-W^^^/J^ (2.56) 

yk ^ jr(k)Q<^Py^ ^ J-W^^^y^, (2.57) 

and W'^^ = 0, which means that the spurious components of one-body operators F^^"^ do not 
play any role in the projected LHA formalism. All we need to do is to keep constant all 
the way along the collective path. Then, we can use the original set of one-body operators 
l^(fc) j on which the LHA equations are projected, without modifying Eqs. ( p.49|) to ( p. 531) . 

In the classical theory of the P+Q model discussed in the next section the constants of 
motion are coordinates g^, and not momenta pj. Therefore this is exactly the case where the 
mixing of spurious components does not play any role. For further details, see the discussion 
around Eq. ( p.21| ). 
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III. COLLECTIVE COORDINATES FOR THE P+Q MODEL 



In this section, we investigate the structure of the self-consistent collective coordinates 
(cranking operators) at the Hartree-Bogoliubov (HB) state for the pairing-plus-quadrupole 
(P+Q) Hamiltonian, as a first step towards the large amplitude collective motion in heavy 
nuclei. We also try to find a set of one-body operators which can approximate the RPA 
eigenvectors. 

A. The quasiparticle RPA for the P-|-Q model 

We apply the symplectic version of LHA with the extended adiabatic approximation to 
the P+Q Hamiltonian. At the HB state, as shown in Sec. [11 (J| , the force condition can be 
always satisfied and the local RPA equation is equivalent to the quasiparticle RPA for the 
constrained Hamiltonian 

H' = H~ Yl ^rNr, (3.1) 



T=n,p 



H = Y: ^AcU - E ^ {PlPr + PrPf) - f E Q\kQ2K (3.2) 
k T=n,p ^ ^ K=-2 

= E (^kClCk - ^ E t^aRaRa + ^ E t^irSaSa, (3.3) 

where are spherical single-particle energies and A^,- = I]fceT c^Cfc are the number operators 
for neutrons (r = n) and protons (r = p). The operators R„ and 5*0- are the hermitian 
and anti-hermitian components, respectively, of the pairing operators, = J2k&T,k>o(^WkJ 
and the dimensionless quadrupole operators, Q2K = ^0 ^ Sfc«(^k^^2_ft:|0clQ, where 60 = 
{h/mujQy^'^ is the harmonic oscillator length. The Hamiltonian contains five operators of 
-Rcr-type and four of S'^-type. Together with the corresponding coupling constants n^, they 
are given by 

Rcr = {P+ ^)ni (-Pf ^)p) ^20^ ^21^ ^22^ 

= (PW)„, (P^))„ Qit\ Qi-,\ (3.4) 



where 



(Pi+^). = -^{Pr ± Pt), for T = n,p, 

= ^(Q2X±Q2-i^), forir = 0,l,2. (3.5) 

The signs (±) indicate the signature quantum number, e'^'^'^^'O^'^^e'^'"'^'' = iO^^-*. Following 
the standard formulation of the model, we will neglect the Fock terms, the contributions of 
the pairing force to the Hartree potential and those of the quadrupole force to the pairing 
potential |p|-P^. After minimising the HB total energy and diagonalising the HB matrix. 
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we obtain the HB ground state |\I/o) and the quasiparticle energies Ek. Using |\I/o) as the 
reference state, the classical Hamiltonian is calculated as 



2 

k cr \i>j ij 



\'>j ij 



Applying the Holstein-Primakoff mapping (2^) and taking into account that the R'{ij) are 



symmetric and S'{ij) are anti-symmetric with respect to the indices (ij), one can see that 
(\l/(t)|-Ro-|\E'(t)) ((\l'(t)|S'o-|\l'(t))) are real (pure imaginary), and indeed the Hamiltonian ( p.6|) 
is real. Differentiating Ti' with respect to vr and ^, the mass and curvature parameters at 
the HB state, defined by Eqs. ( |2.45|) and ( |2.46|) , are given by 



T=n,p 



EJ^p-2Y,K,SM)SM. (3.7) 



V,o.p{i = vr = 0) = - >^rf, 



Nr 

al3 



T=n,p 



= EJ^p - 2 Y ^aRa{a)R^{(3). (3.8) 
(J 

Here, again the Greek index indicates a pair of 2qp index, so that R{a) are the 2qp matrix 
elements R{a) = R{ij) and E^ are the 2qp energies E^ = Ei + Ej. We see that the terms 
R'{ij)pji and S'(ij)pji in the Hamiltonian do not contribute to the mass and curvature 
parameters, which is consistent with the fact that the terms proportional to (a^a) in one- 
body operators are neglected in the quasiparticle RPA. 

Solving the RPA is equivalent to diagonalisation of the matrix B'^'^V^p. For separable 
forces such as the P+Q model, the problem of diagonalisation can be reduced to a root search 
in a multi-dimensional dispersion equation, which facilitates the numerical calculations for 
heavy nuclei [2^]. If the ground state has an axial symmetry, the RPA matrix is block- 
diagonal and can be divided according to K quantum numbers (eigenvalues of Jz). However, 
as we will show below, some nuclei have triaxial HB ground states for which K is no longer a 
good quantum number. Thus, we use the signature quantum number r, which is good even 
for triaxial cases, and the RPA dispersion equation can be divided into positive (r = +1) 
and negative (r = —1) signature sectors. The construction of the response functions and 
the coupled dispersion equations can be found, for example, in Ref. P3 . 



B. Numerical results 



We follow the second and third of the series of papers by Baranger and Kumar |[TT|J12 



for details of the P+Q model. The model space consists of two major shells of A^^osc = 5,6 
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and 4, 5 for neutrons and protons, respectively. Different harmonic oscillator frequencies are 



used for neutrons and protons in order to make the root mean square radii the same |TT 



The single-particle energies are taken from Table 1 in Ref. |I2|, where they are listed in 
units of the harmonic oscillator frequency hujQ = 41.2^4"^/^ MeV. The pairing force strengths 
are Gn = 11A^^ MeV and Gp = 27 MeV for neutrons and protons, respectively. The 
quadrupole force strength is x = TOA"^'^ MeV. The effective charge in the E2 operator 
is taken as e„ = 1.5ZA~^ and Cp = 1 + e„ for neutrons and protons, respectively. Since 



matrix elements of the quadrupole force in upper shells are known to be too strong |]TT 



we 

modify the quadrupole operators by multiplying all quadrupole matrix elements by a factor 
C = (iVL + f)/(^osc + f), where Nose is the oscillator quantum number of the shell under 
consideration and A'^l is that of the lower major shell. Furthermore, due to the difference 
of harmonic oscillator frequencies for neutrons and protons, the operators are multiplied by 
factors = {2N/Ay^^ and = (2Z/A)'^^^ for neutrons and protons, respectively. Thus, 



the dimensionless quadrupole operators in the interaction ( |3.4| ) are actually defined by 

Q2K = {Q2K)n + {Q2k)p, 
iQ2K)n = alChf{r^Y2K)n, {Q2K)r> = alCh^\r^Y2K)p. (3.9) 

These operators will be referred to simply as the "quadrupole operators" in the following 
discussion. 

We have performed calculations for even-even nuclei in the rare-earth region (A^ = 82 ~ 
126, Z = 50 ~ 82) and have investigated the structure of RPA eigenvectors. Since (3 and 7 



vibrations in these nuclei were studied by Baranger and Kumar in Ref. |]12[, we can compare 
our microscopic RPA results with their semi-microscopic ones. 

In Table |, results of the HB and the RPA calculations are presented. The equilibrium 
deformations 7, A„, Ap) at HB ground states are found to agree with Table 2 in Ref. 



|l^ , except ^^''Pt which has a small triaxiality (7 = 52.6°) in our result but has a oblate 



shape (7 = 60°) in theirs. The deformation parameters, (3 and 7, of the HB states are 
defined by [|ll],|ll, 

x(*o|Q5^ 1^0) = ^u;o/5cos7, x(^oIQS^I^o) = huJoP^ini- (3.10) 

We have solved the RPA in the positive signature sector where the spurious modes, the 
pairing and spatial rotation, are coordinates (not momenta) and are characterised by zero 
eigenmodes of the mass matrix, S^^/jj = 0. For spherical nuclei {(3 = 0), the lowest excited 
states with multipolarity A = 2 are selected as quadrupole vibrations and presented in 
columns 7 and 8. For prolate (/5 > 0,7 = 0) and oblate (/5 < 0,7 = 0) deformed cases, 
the first excited states (except zero spurious modes) with K = and K = 2 are taken as 
j3 and 7 vibrations, respectively. However, for some cases, the lowest excited states appear 
to be non-collective and we have selected higher-energy states. These states are indicated 
with * in the table. For a triaxial ground state, the first excited states are selected and 
the excitation energies are shown in column 7 while the transition amplitudes are given in 
columns 8 and 10. The E2 amplitudes M{E2)k are intrinsic values and calculated as 

M{E2)k = I (nW |e„C(r^>^2i.)^^ + e.ar'V^KY,-'^ \0) \ , (3.11) 

where |0) and \n^^^) are the RPA ground and excited states with signature r = +1, re- 
spectively. For spherical nuclei, states with = 0, 1 and 2 are degenerate in energy and 
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we select the K = one for For triaxial shapes, the state 1^*^+)) has non-zero E2 

amphtudes for both K = and K = 2 which are shown in columns 8 and 10 respectively. 

First we compare the RPA frequencies with the harmonic formula given by Baranger and 
Kumar, Eqs. (13), (14) and (15) in Ref. [T^. Unlike our results, this formula contains a core 



contribution to the collective mass (which is denoted by B'^). Since our RPA calculations 
are carried out in the same model space, we should compare the RPA frequencies with those 
calculated by setting = 0. These are shown in Fig. |l]. Open symbols are the results of the 
harmonic formula while solid ones are the RPA frequencies. For spherical and triaxial nuclei, 
the same RPA frequencies are displayed in the figure for f3 and 7 vibrations. The figure shows 
that the frequencies obtained by the formula are larger than the RPA frequencies. For 7 
vibrations, the isotope dependence is well reproduced, even though the absolute values are 
mostly overestimated by up to a factor two. For j3 vibrations, the harmonic formula breaks 
down, especially in well-deformed nuclei. It fails to reproduce both the isotope dependence 
and absolute values. 

Let us try to understand the origin of the harmonic formula. It can be derived by 
taking the small-amplitude limit of the adiabatic time-dependent theory of collective motion 
developed by Baranger and Kumar [|l3l , p!4| which assumes that the collective coordinates are 
the expectation values of the mass quadrupole operators (more precisely the deformation 
parameters of the Nilsson potential). In contrast with the formula, the RPA normal- mode 
coordinates are self-consistently determined by diagonalising the RPA matrix. Therefore, 
the difference between open and closed symbols seen in the figure can be attributed primarily 
to the difference in the collective coordinates. 

Although it is not the purpose of this paper to reproduce the experimental data, we also 
show observed excitation energies of f3 and 7 bands for axially deformed nuclei (Z = 60 ~ 74) 
in Fig. H p8[ . Agreement with the experimental data are significantly improved by choosing 
the collective coordinates properly. One can see that the isotope dependence is roughly 
reproduced in this model. However, we still overestimate the energies by typically about 
500 keV, which is at least partly due to the neglect of anharmonic effects, and may point at 
the limitations of the simple model. 

In Fig. 1^, we show E2 transition amplitudes for the RPA states. Evidently the E2 
amplitudes are smaller for jS vibrations than 7 vibrations. This is partly because the f3 
vibrations are less collective, but also because the collective coordinates for (3 vibrations are 
more complex than those for 7 motion (see below). 

Now let us analyse the structure of the RPA normal-modes by using the projection 
techniques described in Sec. p D| , using various sets of one-body operators. First we utilise 
elementary multipole operators {F^'^^} on which the RPA matrix is projected. We adopt 
monopole operators (rank-0), quadrupole operators with radial dependence r° and (rank- 
2, K = 0,2), those with spin dependence (rank-2, K = 0,2), and hexadecapole operators 
(rank-4, K = 0,2): 

F^^^ = (P+)., (P_)„ iQ20)r, (g22)r, (r^ao)., {r'Y22)r, 

{[r'Y, X s]gLo)., {[r'Y, x sf^l,)^, {[r'V, x s]i^Lo)., {[r'Y^ x s]^^^,)^, 
(r"no)r, (r'>42)r, r = n,p, (3.12) 



where all operators have the positive signature r = +1, as in Eq. ( |3.5| ), but the superscript 
(+) denoting the signature quantum number is suppressed from here on. All matrix elements 
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of rank-2 operators are multiplied by factors ( and a^, just as for the quadrupole operators 
( p.9|) . For /3 vibrations, there are fourteen relevant operators oi K = while there are ten 
operators of K = 2 for 7. For a spherical case, the number reduces to eight (only rank-2 
operators). For a triaxial case, all 24 operators can mix together. 

The calculated frequencies hCl of the projected RPA equation ( p.49| ) are shown in Fig. |^ 



as open symbols with dashed lines and are compared with the real RPA frequencies (solid 
symbols). As is clearly seen, the projection onto the operators ( ^.12[ ) fails to reproduce 
the excitation energies both for (5 and 7 vibrations. Not only the absolute values but also 
the isotope dependence turns out to be incorrect for all nuclei (except for ^°^Pb). The 
result indicates that it is very difficult to obtain sensible approximation by using elementary 
one-body operators. This is mainly due to the fact that the projected RPA eigenvectors 
have unrealistically large amplitudes for high-lying 2qp components. In Ref. 0, we have 
demonstrated for a few Sm isotopes that this can be remedied by introducing a cut-off 
energy for the 2qp matrix elements. For ^°^Pb, which is spherical and has A„ = Ap = 0, the 
projection accidentally works well because the set of available Iplh modes is very limited, 
'^(«ii/2(^i3/2)"^), (5-9/2 (« 13/2)""^), 7r(/i9/2(/iii/2)~"^), T^if 7/2(^1/2)'^) , all of which have almost 
the same particle-hole excitation energies. 

In order to check the effect of suppressing the high-energy components, we also perform 
the projected RPA calculation with the state-dependent operators whose 2qp matrix ele- 
ments are weighted with a factor (ii^2qp)^^- This means that we employ a set of hermitian 
one-body operators {F^'^^} defined by 

^ = E^(«^«^)" + h.c., (3.13) 

where a indicates the 2qp index and F{a) = (q;|F|0). This suppression factor (i?2qp)~^ can 
be derived analytically as follows if the residual Hamiltonian consists only of a single-mode 
separable force H = —^kRR and the RPA frequency is much smaller than 2qp energies. In 
the single-mode case, from Eqs. (|3.7| ) and (|3^ ), the RPA equation becomes 



{{E^y - (hny) /,„ = 2KR{a)Y,Ef,Ri/3)lf,. (3.14) 
Thus, we can analytically determine the RPA eigenvectors fa as 

In the limit of HQ ^ E^, this leads to the quoted suppression factor (£'2qp)~^- 
In the calculation, we adopt only the operators in the separable interactions: 

F^'^) = (P+),, (p_)^^ (g^o)., ((522)., T = n,p. (3.16) 

This results in a six dimensional projected RPA matrix for (3 vibrations and in two di- 
mensional one for 7 vibrations and for spherical cases. Even for triaxial cases, the matrix 
dimension is thus only eight. The results are shown in Fig. ^ by open symbols with solid 
lines. The results are dramatically improved from the previous ones and the obtained fre- 
quencies are almost comparable to those of RPA for both (3 and 7 modes of all nuclei. We 
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can also see the significant improvement over results of the harmonic formula, especially for 
the (5 vibrations (Fig. p. This is because we incorporate the monopole pairing operators 
in Eq. ( p.l6| ) in addition to the quadrupole ones. It is worth noting that the frequency of 
the projected RPA gives an upper limit to the real RPA frequency because it is obtained 
by solving the RPA in a restricted space. On the other hand, this is not necessarily true for 



the harmonic formula of Ref. [0 



Let us examine the eigenvectors of projected RPA equation in more detail. For spherical 
nuclei, the RPA eigenvector has a good isoscalar character and can be well approximated by 

P =Cn{Q2)n + Cp{Q2)p, (3.17) 

with c„ ~ Cp. Here the tilde indicates that the matrix elements include a suppression factor 
as in Eq. ( |3.13| ). A violation of this isoscalar character is seen for ^^^Ba, ^^"Ce, ^^^Nd and 
Pb nuclei in which either the neutrons or protons pairing gap vanishes. The component 
with finite pairing gap is found to be larger than the one with zero gap, so that Cn > Cp 
for Pb and c„ < Cp for ^^^Ba, ^^°Ce and ^'^^Nd. For deformed nuclei, where the collectivity 
of the vibrational states is smaller than for spherical nuclei and the pairing modes can mix 
with the quadrupole ones, the situation is more complex. In Tables |T| and fTT], taking Dy 
and Os isotopes as examples, we show the coefficients of the approximate eigenvectors, / = 
Y,k CkF'^^\ with respect to the set of operators in Eq. ( |3.16|) . The 7 vibrations approximately 
have isoscalar character and the collective coordinate can be approximated by 

P = Cn{Q22)n + Cp{Q22)p, (3.18) 

with Cn ~ Cp. For the (3 vibration, no such simple form can be found, since the properties 
of RPA eigenmodes change from one nucleus to the other. The (3 vibrations for i56-i62-|-)y 
possess a quadrupole nature, Cn{Q 20) n+Cp{Q 20) p, with Cn > Cp, although there are substantial 
contributions from the monopole pairing modes. For ^^^Dy, there is a dominant pairing 
component {P+)n- For Os isotopes, i82~i86Qg j^g^yg dominant proton pairing {P+)p, while 
^^^Os has a quadrupole nature (c^ < Cp). For the triaxial nuclei iss.igoQg^ although the 
K = and K = 2 components are mixed, the vibrational state is of isoscalar quadrupole 
nature with Cn ~ Cp. The changes observed in the RPA normal-mode coordinates for the /3 
vibrations are difficult to be systematically summarised and turn out to be very sensitive to 
details of the underlying shell structure. 

In order to confirm the importance of the monopole components for jS vibrations, we 
have performed a calculation using only the quadrupole operators weighted with the factor 
-^2qp) {Q2o)t and (Q22)t (t = n,p). For spherical nuclei, there is no change from the previous 
results since the monopole and quadrupole components are exactly decoupled. For deformed 
nuclei only, the results are shown in Fig. ^ by open symbols with dotted lines. Differences 
between solid and dotted lines come from the contributions of monopole pairing operators. 
The figure clearly shows that the normal-mode coordinate for (3 vibrations has a significant 
amount of monopole components. 

We now wish to quantify the quality of the approximation. One measure is the difference 
between the real RPA frequencies hfl and the projected RPA hCl, but as usual measures 
based on the eigenvectors are more sensitive. We define two criteria to directly check the 
closeness between the eigenvectors / and /. Choosing the normalisation (of / and /) 
f,aB"^^f^l3 = J^aB^^J^p = 1, those are calculated as 
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Sb = laB'^^ilp - f,p) = 1 - laB'^^U (3.19) 
6i = l- , ^^'^^ , (3.20) 

v/(/, /)(/,/) 

where (/, /') = Y.a f,af'a- We have < 5 < 1 with 5 = corresponding to the exact 
projection and 5 = 1 to the case where / is orthogonal to /. The values of 6b and 6i 
are listed in Tables |1| and |T|. From the values of 6b, we see that the projection works 
reasonably well for Dy and Os isotopes. The worst case is the (3 vibration of ^^^Os for which 
6b ~ 0.3. Roughly speaking, the one-body operators ( |3.16D possess 70% of overlap with the 
real eigenvectors for ^^^Os {6b ~ 0.3), and more than 90% for the others {6b ^ 0.1). 

The values of 6i for f3 vibrations are often much larger than the corresponding 6b (see 
i56,i58£)y and triaxial nuclei ^^^'^^°0s). This indicates that the / for (3 contains some spurious 
components. In the present calculations, the pairing and spatial rotations are the obvious 
spurious motions which are zero eigenmodes of the mass matrix, 

5"'3/^;5 = 0, for f = N^,Np,J,. (3.21) 

The measure 6b (|3.19|) is insensitive to admixture of these spurious components in /. On 
the other hand, the value of 6i will increase by mixing in spurious components. As we have 
discussed in the last part of Sec. |1ID|, the spurious modes do not affect the projected LHA 



equations. Therefore, the large values of 6i should not cause serious problems in application 
of the projected LHA to the large amplitude collective motion. 

In Fig. ^, the values of 6b are shown for the projected RPA using the set of elementary 
operators (|3.12|) (represented by dashed lines) and the set of state-dependent operators 
containing the factor .E^p (|3.16|) (solid lines). For 7 vibrations, the suppression factor E2qp 
can improve the quality of projection by an order of magnitude and we have 6b ^ 0.1 for 
all nuclei. For (3 vibrations, there are some nuclei in the range Z = 68 ~ 80 for which the 
operators in Eq. ( |3.16| ) are not enough to produce the good approximation. This confirms 
the complexity of the normal-mode coordinates for /?. 

We also see a strong isotope dependence in Fig. ^a. For instance, the Er isotopes show a 
staggering behaviour in 6b- This behaviour is found to be closely related to the collectivity 
of the states. Since the (3 vibrations contains the pairing collectivity in addition to the 
quadrupole one, the E2 amplitudes are not necessarily the best indicator of the collectivity 
of the states. We show in Fig. |^ a following quantity for the RPA eigenmodes of f3 vibrations 
in Er nuclei |E6||: 



JMLfT = 3 — ■ — ■ r, n ■, (no summation witli respect to u). (3.22) 

j:,f',9UmM + {m-'9',9',) 

If n 2qp components equally contribute to the RPA mode, we have A^efr = n. Thus, Eq. 
( p.22| ) gives an effective number of 2qp excitations of which the RPA eigenmode consists. 
The figure shows that the collectivity of f3 vibrations in Er isotopes is weak (compared to the 
7 vibrations for which N^s ^ 8 in most cases), and A^cfr shows a staggering pattern correlated 
with that for 6b (while the E2 amplitudes in Fig. ^ show no sign of staggering behaviour). 
For there are dominant 2qp components in the RPA eigenmodes, which decrease 

Nes and increase 6b- For those cases where 6b ^ 0.3 even with state-dependent operators, 
the (3 vibrations are not collective at all, A^^cfT ^ 3. 
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IV. CONCLUSION AND DISCUSSION 



We have examined different version of LHA techniques to see whether the conservation 
laws in the quantum level are satisfied in the classical LHA. It turns out that the restriction 
to point transformations must be lifted in order that spurious motion leads to the exact 
zero modes. The symplectic version of the LHA, which includes an extended adiabatic 
approximation, can guarantee the exact separation of zero-energy spurious modes. 

The projected LHA can be used to truncate the local RPA calculations. We have shown 
that there is no problem with spurious modes for the projection method, at least when the 
spurious modes are zero eigenmodes of the mass matrix. This is exactly the case for the 
pairing and spatial rotation in the P+Q model for which we have investigated the struc- 
ture of the collective coordinates for quadrupole vibrations and examined the possibility of 
expressing the self-consistent cranking operator in terms of a limited number of one-body 
operators. It seems very difficult to approximate the normal-mode vectors with only elemen- 
tary one-body operators. The difficulty disappears, however, when we use a small number 
of state- dependent one-body operators. At HB equilibrium, the collective coordinates for 
quadrupole vibrations in spherical nuclei and for 7 vibrations {K = 2) in axially deformed 
nuclei can be roughly approximated as / ~ Q2K, where Q2K is a state-dependent "mass" 
quadrupole operator defined by Eq. ( ^.9| ) weighted by a factor i?^p for each 2qp matrix 
element ( p.l3| ). However, for (3 vibrations in axial deformed nuclei and for all excitations 
in triaxial nuclei, the monopole pairing components are as important as the quadrupole 
ones. It is worth emphasising that the collective coordinates are very different from the 
usual (state-independent) mass quadrupole operator even for 7 and spherical cases. This 
shows the importance of a self-consistent determination of the collective coordinates for 
large amplitude collective motion, because the coordinates are now found to have a strong 
state-dependence as well. The structure of the self-consistent cranking operators is clearly 
changing when we move from spherical to axially deformed and triaxial nuclei. For the study 
of large amplitude collective motion in heavy nuclei for which the diagonalisation of the RPA 
matrix becomes too time-consuming, the results of this paper gives a suggestion for a choice 
of a state-dependent basis of operators. The choice of a limited set of (state-dependent) 
basis operators provides a practical way to solve the LHA through the projection. With the 
use of self-consistent cranking operators, the LHA should provide a significant improvement 
over a conventional CHFB calculation based on fixed cranking operators. 

Despite all the above qualities, the symplectic version of LHA has a difficulty, namely 
the determination of /^^^, the second derivatives of the final coordinates with respect to 
the original coordinates For the extended adiabatic transformation ( |2.35|) , we also need 
to calculate /(i)^'"^, the terms of second order in tt. At HB states, we need these quantities 
for spurious coordinates /i = / only, /j^^ and /(i)^"^, which can be calculated because we 
know a priori the one-body operators corresponding to the spurious motions. We have done 
this for P+Q model. Away from equilibrium, however, in addition to these = /, we need 
the second derivatives of the collective coordinate, /i = 1. This is not trivial because the 
local RPA equation determines only the first derivatives, f]^ and g^. Unfortunately, so far 
we do not have a general method to determine /]^^ and /(^)^"'^. The projection method may 
actually provide us with a partial answer [^. If we project the LHA equation onto a set of 
one-body operators {F'^'^^}, the equation provides us with a local representation of collective 
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coordinate in terms of a linear combination of the F^'^^^s. Using the exphcit form 

pik) = J2F^'\zj){ala] + a,a,)+J2F'^'\zj)ala„ (4.1) 

i>j ij 

we find 

/i = 1^ = ^ ^^-^-^ = ^EckF^'H^)^ (4-2) 

where F^''\a) = F^^\ij) for 2qp index a = (ij) and we assume that the operators F^''^ are 
hermitian. Here we use a local coordinate system, in which the point under consideration is 
the origin of the phase space, and the relation (2/7) to obtain Eq. (4^). Since we know the 
explicit form of the operators we are now able to calculate the second derivatives f\g 
and second-order terms in momentum 

/i. - ^ - Ec^^S = ^EckG^'\a(3l (4.3) 

f^'^""' - TTir- ^T.ck^p^ = 4Ec,G(^) («/?), (4.4) 

where we neglect the derivative of c^, for simplicity. G'^^^alS) are related to F'^''\ij), by 

G^''\aP) = F'^''\ij)6M, fora = (zfc),/3=(j7). (4.5) 

Although, following this procedure, it is possible to calculate /^^^ and /(i)^"^, we have to 
note that, unlike the determination of f \ by means of the RPA, the determination of these 
quantities through Eqs. ( [4.3|) and ( |4.4| ) does not depend on the dynamics of the system. In 
this approach, part of the dynamics is hidden in the coordinate dependence of the coefficients 
Cfc, which we have ignored for this simple discussion. 

In this paper, we have mainly discussed the structure of the self-consistent cranking 
operators for the P+Q model only at HB minimum points. The structure at points away 
from the minimum and the further application of the LHA technique will be the subject of 
a future publication. 
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TABLES 



TABLE I. HB equilibrium shapes (/3,7), pairing gaps A^, r = n,p (in units of MeV), RPA 
frequencies (in MeV) and E2 transition amplitudes M{E2)k (in efm^). The RPA results for 

spherical nuclei are shown in columns 7 and 8 (columns 9 and 10 are left blank), while for triaxial 
ones the RPA results are in columns 7, 9 and 10 (column 8 is blank). Entries labelled with asterisks 
indicate that they do not correspond to the lowest excited states. See the main text for further 
dclaik. 





Z 


N 


(3 


7 


An 


P 




M(E2)o 




M(E2)2 


Ba 


56 


82 


0.0 


0.0 


0.0 


1.557 


2.87 


30.88 










84 


0.0 


0.0 


0.913 


1.528 


2.24 


35.21 






Ce 


58 


82 


0.0 


0.0 


0.0 


1.686 


2.77 


37.14 










84 


0.0 


0.0 


0.895 


1.654 


2.06 


41.43 










86 


0.0 


0.0 


1.198 


1.624 


1.34 


54.11 






Nd 


60 


82 


0.0 


0.0 


0.0 


1.766 


2.69 


42.58 










84 


0.0 


0.0 


0.878 


1.734 


1.91 


46.91 










86 


0.0 


0.0 


1.176 


1.703 


1.10 


64.37 










88 


0.184 


0.0 


0.964 


1.386 


1.57 


27.51 


1.68 


23.92 






90 


0.236 


0.0 


0.981 


1.229 


1.59 


27.58 


1.75 


24.52 


Sm 


62 


84 


0.0 


0.0 


0.861 


1.778 


1.78 


51.66 










86 


0.0 


0.0 


1.155 


1.746 


0.88 


75.68 










88 


0.198 


0.0 


0.910 


1.394 


1.57 


25.37 


1.63 


24.09 






90 


0.247 


0.0 


0.944 


1.228 


1.42 


28.97 


1.69 


23.83 






92 


0.312 


0.0 


0.902 


1.007 


0.66 


34.58 


1.90 


21.71 


Gd 


64 


84 


0.0 


0.0 


0.845 


1.793 


1.68 


55.34 










86 


0.0 


0.0 


1.134 


1.761 


0.71 


87.45 










88 


0.203 


0.0 


0.873 


1.386 


1.55 


23.66 


1.59 


24.61 






90 


0.249 


0.0 


0.917 


1.221 


1.44 


27.33 


1.62 


25.65 






92 


0.299 


0.0 


0.904 


1.056 


0.89 


29.62 


1.69 


26.74 






94 


0.320 


0.0 


0.880 


0.980 


1.15 


22.52 


1.60 


27.64 






96 


0.329 


0.0 


0.840 


0.938 


1.40 


17.82 


1.35 


26.36 


Dy 


66 


90 


0.244 


0.0 


0.897 


1.212 


1.55 


24.78 


1.58 


27.71 






92 


0.285 


0.0 


0.901 


1.069 


1.22 


25.33 


1.52 


29.88 






94 


308 





873 


988 


1 3Q 


20 09 


1 45 


31 75 






96 


0.320 


0.0 


0.827 


0.941 


1.56 


13.31 


1.23 


30.98 






98 


0.327 


0.0 


0.748 


0.904 


1.50 


0.59 


1.13 


29.35 


Er 


68 


94 


0.295 


0.0 


0.867 


0.986 


1.58 


20.12 


1.49 


30.81 






96 


0.310 


0.0 


0.817 


0.927 


1.59 


9.18 


1.28 


29.50 






98 


0.320 


0.0 


0.736 


0.884 


1.47 


1.98 


1.17 


27.45 






100 


0.326 


0.0 


0.680 


0.849 


1.36 


1.03 


1.51 


26.37 






102 


0.330 


0.0 


0.624 


0.819 


1.20 


5.40 


1.72 


26.95 


Yb 


70 


94 


0.278 


0.0 


0.874 


1.016 


1.59 


20.62 


1.56 


27.79 






96 


0.296 


0.0 


0.817 


0.946 


1.58 


11.28 


1.40 


24.94 






98 


0.309 


0.0 


0.734 


0.892 


1.45 


4.82 


1.29 


21.61 






100 


0.316 


0.0 


0.679 


0.852 


*1.66 


*11.33 


1.67 


17.41 






102 


0.320 


0.0 


0.632 


0.820 


1.18 


8.21 


*2.27 


*22.36 
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104 

106 

Hf 72 94 
96 
98 
100 
102 
104 
106 
108 
110 

W 74 98 
100 
102 
104 
106 
108 
110 
112 

Os 76 106 
108 
110 
112 
114 
116 

Pt 78 106 
108 
110 
112 
114 
116 
118 
120 

Hg 80 112 
114 
116 
118 
120 
122 
124 

Pb 82 118 
120 
122 
124 
126 



0.313 


0.0 


0.305 


0.0 


0.252 


0.0 


0.274 


0.0 


0.289 


0.0 


0.299 


0.0 


0.302 


0.0 


0.293 


0.0 


0.282 


0.0 


0.272 


0.0 


0.257 


0.0 


0.264 


0.0 


0.275 


0.0 


0.277 


0.0 


0.270 


0.0 


0.259 


0.0 


0.247 


0.0 


0.228 


0.0 


0.202 


0.0 


0.235 


0.0 


0.222 


0.0 


0.201 


0.0 


0.181 


22.1 


0.171 


46.0 


-0.156 


0.0 


0.205 


0.0 


0.190 


0.0 


0.171 


20.1 


0.160 


52.6 


-0.149 


0.0 


-0.138 


0.0 


-0.126 


0.0 


-0.107 


0.0 


-0.139 


0.0 


-0.131 


0.0 


-0.123 


0.0 


-0.114 


0.0 


-0.089 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



0.739 


0.813 


0.803 


0.811 


0.907 


1.060 


0.842 


1.004 


0.759 


0.963 


0.706 


0.932 


0.678 


0.907 


0.744 


0.893 


0.785 


0.884 


0.824 


0.878 


0.920 


0.883 


0.819 


0.962 


0.768 


0.940 


0.752 


0.918 


0.777 


0.893 


0.796 


0.872 


0.838 


0.856 


0.932 


0.854 


0.988 


0.872 


0.838 


0.785 


0.883 


0.762 


0.966 


0.758 


0.990 


0.789 


0.982 


0.830 


0.928 


0.844 


0.930 


0.754 


0.976 


0.723 


1.046 


0.671 


1.078 


0.674 


1.030 


0.676 


0.944 


0.676 


0.814 


0.678 


0.693 


0.692 


1.106 


0.0 


1.045 


0.0 


0.956 


0.0 


0.830 


0.0 


0.745 


0.349 


0.749 


0.606 


0.510 


0.596 


1.011 


0.0 


0.894 


0.0 


0.734 


0.0 


0.497 


0.0 


0.000 


0.0 



*1.55 


*13.39 


1.42 


18.82 


1.63 


18.96 


1.65 


10.49 


1.47 


8.59 


*1.72 


*15.93 


1.22 


12.52 


*1.69 


*14.97 


1.56 


4.48 


1.66 


0.28 


1.54 


27.21 


1.53 


13.74 


1.47 


13.95 


1.31 


16.48 


1.53 


8.60 


1.59 


0.68 


1.63 


9.29 


1.71 


16.35 


1.70 


21.51 


1.57 


2.00 


1.51 


7.88 


1.47 


12.34 


0.59 


18.96 


0.52 


44.03 


1.64 


17.81 


*1.66 


*17.68 


1.44 


6.59 


0.72 


17.05 


0.21 


72.21 


1.44 


6.93 


1.50 


6.76 


1.49 


18.53 


0.99 


26.14 


1.61 


17.73 


1.69 


27.36 


1.60 


26.48 


1.52 


18.62 


0.43 


32.92 


0.51 


49.51 


0.90 


29.38 


0.80 


36.93 


0.89 


29.92 


1.01 


22.41 


1.12 


14.14 


4.71 


19.86 



2.02 


25.66 


1.68 


27.51 


1.43 


32.01 


1.38 


28.52 


1.33 


23.97 


1.68 


19.10 


*2.32 


*23.19 


1.99 


26.79 


1.60 


29.93 


1.30 


34.17 


1.03 


40.07 


1.28 


33.44 


1.60 


29.98 


1.95 


34.85 


1.73 


36.22 


1.40 


38.14 


1.11 


43.06 


0.79 


52.97 


0.37 


83.11 


1.12 


47.79 


0.85 


54.86 


0.42 


79.52 




63.99 




40.94 


0.61 


53.03 


0.68 


63.59 


0.23 


110.48 




55.93 




55.38 


0.44 


58.97 


0.83 


40.18 


1.20 


30.09 


1.29 


21.97 


0.88 


35.81 


1.00 


30.80 


1.23 


24.73 


1.53 


17.54 


1.37 


13.22 



23 



TABLE II. RPA frequencies HQ, projected RPA frequencies HQ, (both in units of MeV) and 
projected RPA solutions for (3 and 7 vibrations in Dy isotopes. Coefficients for the approximated 
RPA normal-mode coordinates, Ck of f^a = J2k'^kF^''\c(), are listed normalized as J2k kfcP — 1- 



The two possible 


measures 


of the quality of projection, 6b and 5i 


, are shown at the bottom. 




156 


Dy 


158 


Dy 


160 


Dy 


162 


Dy 


164Dy 




P 


7 




7 


(3 


7 


(3 


7 


(3 


7 


nn 


1.55 


1.58 


1.22 


1.52 


1.39 


1.45 


1.56 


1.23 


1.50 


1.13 


hn 


1.77 


1.77 


1.34 


1.65 


1.53 


1.54 


1.85 


1.30 


1.94 


1.20 


iP+)n 


-0.061 


0.0 


-0.071 


0.0 


0.030 


0.0 


0.106 


0.0 


0.988 


0.0 


{P-)n 


-0.494 


0.0 


-0.437 


0.0 


-0.350 


0.0 


-0.312 


0.0 


0.064 


0.0 


(<52o)n 


0.609 


0.0 


0.671 


0.0 


0.745 


0.0 


0.784 


0.0 


0.142 


0.0 


(<522)n 


0.0 


0.730 


0.0 


0.706 


0.0 


0.693 


0.0 


0.716 


0.0 


0.726 


{P+)p 


-0.325 


0.0 


-0.246 


0.0 


-0.225 


0.0 


-0.240 


0.0 


-0.009 


0.0 


{P-)p 


-0.026 


0.0 


-0.024 


0.0 


-0.021 


0.0 


-0.019 


0.0 


-0.001 


0.0 


{Q20)p 


0.525 


0.0 


0.540 


0.0 


0.520 


0.0 


0.467 


0.0 


0.016 


0.0 


{Q22)p 


0.0 


0.684 


0.0 


0.708 


0.0 


0.720 


0.0 


0.698 


0.0 


0.688 


Sb 


0.058 


0.064 


0.010 


0.028 


0.011 


0.008 


0.132 


0.007 


0.137 


0.010 


Si 


0.410 


0.059 


0.406 


0.025 


0.140 


0.004 


0.176 


0.005 


0.115 


0.008 



TABLE III. Similar to Table III but for Os isotopes. Note that the nuclei i88.i90Os have a 
triaxial shape in the ground state. 





1820 


s 


184 


Os 


186 


Os 


1880s 


190Os 


192 


Os 




f3 


7 




7 


P 


7 






/3 


7 


nn 


1.57 


1.12 


1.51 


0.85 


1.47 


0.42 


0.59 


0.52 


1.64 


0.61 


nn 


1.85 


1.15 


1.78 


0.86 


1.73 


0.43 


0.61 


0.53 


2.02 


0.62 


iP+)n 


-0.014 


0.0 


0.265 


0.0 


0.174 


0.0 


-0.014 


-0.125 


0.236 


0.0 


{P-)n 


0.001 


0.0 


-0.081 


0.0 


-0.162 


0.0 


0.208 


0.113 


-0.082 


0.0 


(<520)n 


0.009 


0.0 


-0.186 


0.0 


-0.278 


0.0 


0.178 


0.504 


0.561 


0.0 


(<322)n 


0.0 


0.690 


0.0 


0.699 


0.0 


0.705 


0.658 


0.472 


0.0 


0.705 


iP+)p 


0.999 


0.0 


0.905 


0.0 


0.847 


0.0 


-0.049 


-0.047 


0.381 


0.0 


iP-)p 


0.031 


0.0 


-0.069 


0.0 


-0.101 


0.0 


0.099 


0.064 


-0.117 


0.0 


{Q20)p 


-0.028 


0.0 


-0.256 


0.0 


-0.372 


0.0 


0.184 


0.509 


0.681 


0.0 


{Q22)p 


0.0 


0.724 


0.0 


0.715 


0.0 


0.709 


0.668 


0.479 


0.0 


0.710 


Sb 


0.072 


0.002 


0.067 


0.0005 


0.068 


0.00005 


0.001 


0.001 


0.313 


0.0004 


Si 


0.064 


0.001 


0.087 


0.0003 


0.156 


0.00004 


0.835 


0.657 


0.329 


0.0003 
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FIG. 1. Calculated excitation energies in units of MeV, of P (a) and 7 vibrations (b) as functions 
of mass number A for even-even nuclei in the rare-earth region. The closed circles indicate the 
RPA results while the open symbols are results of the harmonic formula by Baranger and Kumar 
|12|. The frequency of the harmonic formula for ^'^^Pb is 10.9 MeV. For cases in which the HB 
ground states are spherical or triaxial, we cannot distinguish /3 from 7 vibrations and the same 
values have been used for both (a) and (b). Different symbols indicate different isotopes; circles 
for Ba, Dy and Os; squares for Ce, Er and Ft; diamonds for Nd, Yb, Hg; triangles (upwards) for 
Sm, Hf and Pb; triangles (downward) for Gd and W. The experimental values are represented by 
crosses for deformed nuclei 
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FIG. 2. Calculated intrinsic E2 transition amplitudes between the excited RPA and the ground 
state, for /? (a) and 7 vibrations (b) in units of efm^. For spherical nuclei, these values are also the 
E2 amplitudes in laboratory frame. For triaxial cases, a single state has both K = and K = 2 
components of E2 amplitudes which are shown in (a) and (b), respectively. See caption of Fig. |l| 
for symbols and the text for further details. 
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FIG. 3. Calculated energies in units of MeV, by means of both the projected RPA (open 
symbols) and the exact RPA (closed) for P (a) and 7 vibrations (b). The results of projected RPA 
connected with dashed, dotted and solid lines indicate different sets of one-body operators being 
employed for the projection. See the main text for further details. 
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FIG. 4. Quality of projection 5b, Eq. ( p.19 ), for f5 (a) and 7 vibrations (b). Dashed lines and 
solid lines denote different sets of one-body operators being employed for the projection. See the 
main text for further details. 
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FIG. 5. The measure 5b (open squares and scale on the left) and the measure of collectivity 
A'eff, Eq. ( 3.221 ), (closed squares and scale on the right) for f3 vibrations in the Er isotopes. 
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